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Abstract. This review is dedicated to recent results on the 2d parabolic-elliptic Patlak- 
Keller-Segel model, and on its variant in higher dimensions where the diffusion is of critical 
porous medium type. Both of these models have a critical mass M c such that the solutions 
exist globally in time if the mass is less than M c and above which there are solutions which 
blowup in finite time. The main tools, in particular the free energy, and the idea of the 
methods are set out. A number of open questions are also stated. 



1. Biological background 

Chemo-taxis is defined as a move of an organism along a chemical concentration gradient. 
Bacteria can produce this chemo-attractant themselves, creating thus a long-range non-local 
interaction between them. We are interested in a very simplified model of aggregation at 
the scale of cells by chemo-taxis: myxamoebaes or bacterias experience a random walk to 
spread in the space and find food. But in starvation conditions, the dictyostelium discoideum 
emit a chemical signal: cyclic adenosine monophosphate (cAMP). They move towards a 
higher concentration of cAMP. Their behaviour is thus the result of a competition between 
a random walk-based diffusion process and a chemo-taxis-based attraction. It was noticed 
experimentally that if there are enough bacteria they aggregate whereas if they are not enough 
they go on spreading in a chemo-tactically inert environment, e.g. [15] , The typical time scale 
for the spreading of bacteria on a petri dish is around one day, and a few minutes for the 
concentration. This concentration phenomenon is the first step for uni-cellular organisms to 
come together with others and form a multi-cellular organism. It can be seen as a hint on 
how, during the evolution of species, the passage from uni-cellular organisms to more complex 
structure was achieved. It is also a paradigm model for pattern formation of cells for meiose, 
embryo-genesis or angio-genesis. 

In nature the dictyostelium discoideum spread on the soil and then come together to form a 
motile pseudoplasmodium. This slug creeps to a few centimeters below the soil surface where 
it forms a fruiting body with spores and a stalk. The spores are then blown away by the wind 
to colonise a new place. Around 20% of the cells which are in the stalk altruistically sacrifice 
themselves to allow the species to survive. They are an excellent example of social behavior 
with outstanding coordination and sense of sacrifice for the benefit of the species. 
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2. The 2d parabolic-elliptic Patlak-Keller-Segel system 

2.1. The model. The first mathematical attempt to model this aggregation phenomenon is 
often granted to Evelyn Fox Keller and Lee A. Segel in [38] but the same model was earlier 
described by Patlak in [19]. We consider the following simplified version was given by [36 j : 

^ = Ap-V-(pVc) i£R 2 ,t>0, 

-Ac = p x E M 2 , t > 0, W 

^ p(-,t = 0) = po > xeR 2 . 

Here p represents the cell density and c the concentration of chemo-attractant. 

As the solution to the Poisson equation —Ac = p is given up to a harmonic function, we 
choose the one given by c = G * p where 

G{\x\) := -— log\x\ . 

The Patlak-Keller-Segel system ([1]) can thus be written as a non-local parabolic equation: 

^ = Ap - div(pVG * p) in (0, +oo) x E 2 . 
Also note that the mass is conserved 

p(t,x) dx = Po{x) dx =: M . 



S. Childress and J. Percus and V. Nanjundiah conjectured in [261 W7\ that this system 
displays the existence of a critical mass above which the cells aggregate and below which they 
do not., For a complete review of the early literature, the interested reader could beneficially 
consult 1351 134|. For a more recent references see |50| or 



Except when it is clearly indicated, in this article we will assume that the initial condition 
satisfies: 

(1 + |x| 2 )p G L\{R 2 ) and p logp € L 1 ^ 2 ) . 



2.2. Blowup. Consider a smooth solution to the Patlak-Keller-Segel system ([T]) we formally 
compute a virial identity 

d 
~dl 

So that all the solutions with mass bigger than 8ir and finite 2-moment cannot be global in 
time: 



/ \x\ 2 p(t, x) dx = 4M(l - — ) . (2) 

7R2 87T' 



Theorem 1 (Blowup). Let p be a solution to the Patlak-Keller-Segel system (pQ) and [0,T*) 
its maximal interval of existence. If M > 8tt, then 

T ^M(M-8n)Jj X ^ x)dX > 

and p(t, ■) converges, up to extraction of a sub-sequence, as t —> T* to a measure which is not 
in LW). 
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This virial method was already used in the dispersive equations community. However, this 
kind of method is non-constructive and gives no hint on the nature of the blowup. Another 
non-constructive proof of this non global-in-time existence will be given in Remark [JJ and will 
give similar results even in the case when the solutions are of 2-moment infinite. 

The proof that the blowup occurs in an isolated point was initiated in |46t 143] . When the 
solution is radially symmetric, the blowup mechanism is not clear and the different conjectures 
are as follows: There exists a solution to the Patlak-Keller-Segel system ([1]) that yields a 
concentration of an amount of mass 8ir at the origin. Moreover, the asymptotics of this 
solution near the origin is given by 



p{t,x) 



1 



£(*) 2 (i + M 2 ) 2 



and c(t,x) = — 21og(l + y 2 ) 



(3) 



for t — > T~ where y = r/L(t) and as t goes to oo. The different approach lead to a different 
expression for L. The most recent approach is due to [42] where the author makes strong 
analogies between the Patlak-Keller-Segel system and the non-linear Schddinger equation. We 
will come back to these analogies in the concluding remarks in Section P. Lushnikov uses 
an appropriate gauge transformation to a new time dependent variable around the self-similar 
solution to obtain a linearisation operator of self-adjoint form. The study of the spectrum of 
the linearised operator allows him to derive a set of amplitude equations for the coefficients 
of this expansion. He thus obtain 



L(t) = 2\/T-te 



-1-7/2 ^/-log(r-t)/2 



1 + 



[io g (-iog(r-t))r 

y/-2log(T-t) 



(4) 



7 being the Euler constant. In [59], a similar result was obtained using matched asymptotic 
expansion techniques but with exponent 1 in the log- log term in @. Whereas 0T] conjectures 



L{t) ~ VT^eV- 1 ^-')/ 2 (-log(T - 4 ))(-lo^-*))- 1/2 /4 



If the numerics performed in [52] agree on the leading order term y/T — t, C. Sire and P.- 
H. Chavanis predict that 

Lit) ~ y/T^t eT^I- W log(- iog(T-t))/2/2 _ 

Open question: All those estimates are formal. Justification of this result is very challeng- 
ing, as is the numerical track of the next order term in L. 

Concerning the continuation of the solution after blowup, the usual idea is to define a 
sequence of approximate problems containing a small parameter e > which have global- 
in-time solutions and approach the original problem when e goes to 0. The behaviour of 
the solution to the approximate model has to be close to the one of the Patlak-Keller-Segel 
system ([1]), except close to the singularities. The blowup indicates that the approximate 
problem is not valid anymore close to the singularity. In [61|, [60] , using matched asymptotic 
expansions J. Velazquez describes in a rather detailed manner the formation and motion of 
some regions where the mass concentrates. He also proves local existence of the solution as 



1 



ADRIEN BLANCHET 



long as there is no formation of another Dirac mass and no collision of Dirac masses occur. 
The approximate problem is 

^ = Ap-V-($ e (p)VG*p) x G R 2 , t > , 

where <& e (p) '■= s^ 1 ^{ep) where e is a small parameter and $ is an increasing bounded 
function satisfying 

<3?(s) = s — as 2 + • • • as s — > and <3?(s) ~ L as s — > oo , 

where L > is a given number. This model prevents overcrowding as the chemo-tactic 
function & e (p) saturates at the constant value L. 

In [30], J. Dolbeault and C. Schmeiser define measure valued densities to give a sense to 
generalized global-in-time solutions for any mass. This extends the solution concept after 
blow-up. They also show that the choice of a solution concept after blow-up is not unique 
and depends on the type of regularization. Their regularisation is different from the one 
chosen above. The regularised problem they consider is 

^ = Ap- V- {pVG e *p) xeR 2 ,t>0, (5) 
where G e (x) := — log(|x| +e)/2ir. 

Theorem 2 (Generalised solution, [30]). For every T > 0, as e — > 0, a sub-sequence of 
solutions e to ([5]) converges tightly and uniformly in time to a time dependent measure p(t). 
There exists u(t) such that (is, p) is a generalized solution in the distributional sense of 

^+div(j[p,z,]-V P )=0 (6) 
where the convective flux j [p, v\ is supported in the support of v and is given by 



1 



tp(t, x)j[p, u](t, x) dx dt = — — j j u(t, x)V(p(t, x) dx dt 

T 



o JR 2 47r Jo 



T 



1 

4-7T 

for any tp G C£((0, T)xK 2 ) with 



o 



{<p(t, x) - (p(t, y)) K(x - y)p(t, x)p(t, y) dx dy dt 



for x 7^ 
for x = . 



Ixl 2 



If p does not charge points then the additional default measure v vanishes and j[p, 0] = 
pVG* p, so that §6§ is generalisation of the Patlak-Keller-Segel system ([T]). They also obtain 
a strong formulation when the generalised solution is assumed to be the sum of a regular part 
and of Dirac masses: 

Theorem 3 (Strong formulation, |30j). Assume that the generalised solution to ([6]) has the 

form 

p(t, X ) = ^2Mj(t)5(x -Xj(t)) + Prc g {t, x) . (7) 

jeN 
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Then as long as there is no formation of a new Dirac mass and that there is no collision of 
Dirac masses 



dp 



dt 



with 
and 



/ N 
Ap re g - V PregVG * p rcg - ^Vp rcg ^ M 3 ft) i _ ~ J t ) i 2 V ^ G § 

Mi(t) = Mi{t)p Teg (t,Xi(t)) 

i l (t) = VG* Prcg (t,x l (t))-— E M iW uT^2 • 

A similar result was formally obtained in [61], with the last equation replaced by 

Xi (t) = r(M i (t))(-±- £ Mj(t) X J X ^ 2 +VG, Prcg (t,x t (t)) 

where T(M) is a mean value of the derivative <3? is such that < T(M) < 1, r(87r) = 1 and 
r(oo) = and can be explicitly described, see |61[ Equation (3.45)]. Note that this means 
that the limit depends on the choice of the regularisation. 

In [21], for the case of two particles, the authors prove that when t is large enough the 
solution is made of a Dirac peak of mass Ma(t) surrounded by a dilute halo containing 
the remaining mass whose dynamical evolution is described by a Fokker-Planck equation. 
Therefore, they neglect the self-gravity of the halo and prove that the mass of the Dirac peak 
saturates to M algebraically rapidly as 

i M (t) 

Open question: Actually the assumption (|7|) is valid only between two blowup events or 
between two collisions of Dirac masses. We expect that in the end the solution is made of one 
Dirac mass with all the mass but we are still missing such a rigorous theory. 

2.3. Global existence. 

2.3.1. A priori estimates. The first important answer to Nanjundiah's conjecture on global 
existence vs blowup was given by W. Jager and S. Luckhaus in |36| : the natural idea is to 
regularise the Green kernel and to pass to the limit. The solutions have mass M so that the loss 
of compactness can come either from concentration or vanishing. By ([2]), the second moment 
remains bounded so that the main problem is to control the concentration of mass. W. Jager 
and S. Luckhaus tried to obtain a bound on the entropy J R2 plog/5 dx by differentiating it 
and using an integration by parts and the equation for c, we obtain: 

-j- / P log p dx = — 4 f |V-y/p| 2 dx + f p 2 dx . 

dt J R 2 J R 2 J R 2 

Applying the Gagliardo-Nirenberg-Sobolev inequality: 

/ M 4 dx < C GNS / \Vu\ 2 dx [ \u\ 2 dx V u € H 1 (M. 2 ) . (8) 

JR 2 JR 2 Jm 2 
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to u = y/p, we have 

/ plogpdx< [-4 + MCgns] / |V^p| 2 dx 

CIC Jr 2 J R 2 

So that the entropy is non-increasing if M < 4/Cqns ~ 1-862... x 47r < 871". They hence 
obtained global existence in this case together with important propagation of the L p -estimates. 

We can indeed improve this result by using the following free energy: 

•TtkS [/>]:= / plogpdx-i/ pcdx. 
Jm 2 z Jm. 2 

A simple formal calculation shows that for all u S C^ C (M 2 ) with zero mean, 

lim - (J"pks[p + eu] ~ FpKslp]) = [ ^ PKS P (x) u{x) dx 
e->0 e J R 2 dp 



where 

^pks [p] 



(x) := log p(x) — G * p(x) 



bp 

It is then easy to see that the Patlak-Keller-Segel system ([1]) can be rewritten as 

'STp KS [p(t)} 



^(t,x) = div\ 



■ x) 



Sp 

It follows that at least along well-behaved solutions to the Patlak-Keller-Segel system (P), 



(9) 



d f 
-r^PKsW)} = - / P (t, 
dt J R 2 

Or equivalently 



V 



5p 



dx . 



_d_ 



FpKs[p(t,-)]=- I HV(logp-c)| 2 dx. 
Jm. 2 



In particular, along such solutions, 1 1— > ^PKs[p(i)] is monotone non-increasing. 

The gap between the 4/Cgtvs an d 8ir was not filled before [29] when was made the link 
with the logarithmic Hardy-Littlewood-Sobolev: Let / be a non-negative function in L 1 (M 2 ) 
such that / log / and /log(l + |x| 2 ) belong to L^R 2 ). If / R2 f dx = M, then 

/ flogfdx + ^-ff f(x)f(y)log\x-y\dxdy>-C(M), (10) 
Jm. 2 M J ii 2 xi 2 

with C(M) := M(l + log7r — log M). Moreover the minimisers of (|10p are the translations of 
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Using the monotonicity of J-pks [p] an d the Logarithmic Hardy-Littlewood-Sobolev inequal- 
ity it is easy to see that 



p(x) log \x — y\p(y) dx dy 



^pks[p] = ^(/ k/ (x) lo gP (x)dx + ^j^ 
+ n. - ^-Q p(x) \ogp{x) dx 
>~ C(M) + - p(x) log p(x) dx . (11) 

It follows that for solutions p of the Patlak-Keller-Segel system (pQ), 

p( t ,x)log P (M)dx< 8 ^W (12) 

07T — M 

Therefore, for M < 8ir, the entropy stays bounded uniformly in time. This precludes the 
collapse of mass into a point mass for such initial data. 

Coming back to the super-critical mass case, it is worth noticing that for a given p, if we 
set p\(x) = X~ 2 p(X~ 1 x) then 

^pks [pa] = Fpks W -2M\l-j£\ log A . (13) 

So that as a function of A, -7~pks[pa] is bounded from below if M < 8ir, and not bounded from 
below if M > 8ir. 

As an alternative to the regularisation/passing to the limit procedure, another conceited 
but smart way to prove the global existence is to use the gradient flow interpretation in 
the Wasserstein metric. For this purpose we need to introduce a few elements of optimal 
transport. We denote V{X) the set of probability measure on X and V2(X) the subset of 
probability measures with finite second moments. We say that T transports p onto u, and 
denote T#pL = u, if 

C(T(x)) dp(x) = [ C(V) du(y) VC € C° b (R 2 ) . 
We can define the 2- Wasserstein distance, W2 in V(X) x V(Y) by 



ner 



W 2 (/^r=inf// \x-y\ z dU(x,y) , 



XxY 



where II runs over the set of joint probability measures inXxF with first marginal p and 
second v. 

Using this metric, we can see the Patlak-Keller-Segel system JI| as a gradient flow of the 
free energy in the Wasserstein metric: 

Pt = -Vty-FpKsIpW] • 
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In the sense that we can construct a solution using the minimising scheme, often known as 
the Jordan-Kinderlehrer-Otto (JKO) scheme: given a time step r, we define the solution by: 



p T ^ G argmm peJC 



+ -^PKS [p] 



2r 

where K. := {p : f R2 P = M, L 2 p{%) ^ogp{x) dx < oo and J* R2 \x\ 2 p(x) dx < 00} . 

For the analogy, note that if the metric was Euclidean, the Euler-Lagrange equation asso- 
ciated to 



k -1-1 

p T 6 argmm 



2 +^PK S [p] 



2r 



(14) 



would be 



fe+i _ „fe 
^ ^+J-PKS^ +1 ]=0, 

r 

which is nothing but the implicit Euler scheme associated to 

Pt = -VJ=PKs\p{t)] . 

At this point it is convenient to emphasise that the functional J-pks is n °t convex, so even 
the existence of a minimiser is not clear. When the functional is convex, or even displacement 
convex, general results from \62\ Q] can be applied. In concrete terms, a functional Q is said 
to be displacement convex when the following is true: for any two densities po and p\ of the 
same mass M, let (p be such that Vip#po = p±. For < t < 1, define 

cp t (x) = (1 - + and = ^<Pt#Po ■ 

The displacement interpolation between po and p\ is the path of densities 1 1— >■ p tj < t < 1. 
Let 7 be any real number. To say that Q is ^-displacement convex means that for all such 
mass densities po and pi, and all < t < 1, 

{i-t)g(p ) + tg(p l )-g( Pt )> 1 t(i-t)wl(po,p l ) . 

Q is simply displacement convex if this is true for 7 = 0, and Q is uniformly displacement 
convex if this is true for some 7 > 0. 

We interpolate between the terms of the sequence {p^-}keN to produce a function from 
[0,oo) to L 1 (M 2 ): For each positive integer k, let \7ip k be the optimal transportation plan 
with Vtp k #p k = Pt' 1 - Then for (k — l)r < t < kr we define 

M«)-( '~ ( * T ~ 1)T M + *^v)#rf. 

Theorem 4 (Convergence of the scheme as r — )• 0, [8]). J/M < 8-7T then the family (p T )r>o 
admits a sub-sequence converging weakly in L 1 (M 2 ) to a weak solution to the Patlak-Keller- 
Segel system. 

In this proof the perturbation of the minimiser has to be done in the Wasserstein setting: 
Let C be a smooth vector field with compact support, we introduce T e := id + eQ. We define 
~pl the push-forward perturbation of by T £ : 

p- = T £ #p n T +1 . 
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Standard computations give 

f C{x f-^ n (-) p n + l {x) dx 



, i t lew -CM1 •('-») „^ {y)iy 

4vry R2 \x-y\ z 



x) dx 



which is the weak form of the Euler-Lagrange equation: 



id Vy" +1 = ^ + p n+l Vc n+l (15) 

r 

Using the Taylor's expansion £(x) — £ [V<£>"(x)] = [x — Vy? n (x)] • V£(x) + O |x — Vip n (x)\ 2 
we obtain for allt2 > h > 0, 

/ C(^)[Pr(*2,^) -p T (ti,a;)] dx = / / AC(x)p T (s,x) dxds + O^ 1 / 2 ) 
JR 2 Jti JR 2 



1 [ t2 [[ f , ( v (x - • (VC(x) - vcfa)) , , 

— / // p T (s,x)pr{s,y) : j= dydx. (16) 

4vr7 tl JJm.2 xR 2 \x-y< z 



To pass to the limit, the scheme provides some a priori bounds: Taking as a test function 
in (|14p we have: 

^PKSK +1 ] + ^Wl{p n r,p n r +1 ) < ^PKsK] • (17) 
As a consequence we obtain an energy estimate 

sup Jpks[p?] < Jtks [Pt], 
which together with (jlip forbids the concentration, and a totaZ square estimate 

E W 2(P?, < ^pks[p?] - inf -FpksK] , 

which rules out vanishing. These two estimates allow to pass to the limit in r in (|16|) . to 
obtain: 

C(x)[p(t 2 ,x) - p(ti,x)\ dx = / AC(x) p(s,x) dx ds 



4vrJ tl JJ R 2 xR 2 \x-y 



Which is the definition of a weak solution. 

We can actually obtain: 

Theorem 5 (Existence of solution is the subcritical case, [13]). If M < 8ir, then the Patlak- 
Keller-Segel system has a global weak non-negative solution p with initial data po such 
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that 

(l + |x| 2 + |logp|) /0 GL£(]R+,L 1 (M 2 )) 



/ / p\ V log p — Vc| 2 dx dt < oo , 
Jo Jm 2 

p(t,x)dx= / \x\ 2 p (x) dx + AM (l )t 

Jr 2 ^ 



for t > 0. Moreover p € L?° ((e, oo), L' P (R 2 )) /or any p G (l,oo) and any e > 0, and the 
following inequality holds for any t > 0: 

J r PKs\fi(-,t)]+ / p|V (logp- c)| 2 dx ds < -FpksI/Oo] • 
J JR 2 

Similar results were first proved in [45] for radially symmetric solutions in a bounded 
domain with Neumann boundary conditions. 

This notion of free energy solution allows to study the large time behaviour, intermediate 
asymptotics and convergence to asymptotically self-similar profiles: let (woo? 

t>oo) be the unique 

solution to the Gelfand equation 

g«oo-M 2 /2 

Uoo = M- _ 2 = -Avoo , with v O0 = G*u OQ . (18) 

J R2 e °° dx 

Using the comparison principle in the radial variable, it has been proven in [B] that the 
radial, non-negative smooth solution to this problem is unique. In the original variables, the 
self-similar solutions of ([TJ take the expression: 

1 



poo(t,x) := Y^ t u °° (log(v / l + 2t),x/ v / l + 2i) , 



Coo (t,x) := Voo (log(v / l + 2t),x/ v / l + 2t) . 
Theorem 6 (Large time behaviour, [13]). Under the assumptions in Theorem^ 

lim \\p(-,t) - Poo(-,*)||li(r2) = and lim ||Vc(-, t) - Vc^-, i)||z,2( R 2) = . 

The proof follows the usual entropy/entropy production method in self-similar variables: 
We define the rescaled functions u and v by: 

p{t > x) = iW) u {wy T{t) ) and c{t > x) = v {wr T{t) 

with R(t) = yl + 2t and r(t) = log R(t). The rescaled system is 
( du 

— = Au-V ■ (u(x + Vv)) xeR 2 ,t>0, 

v = G*u x£R 2 , t>0, ( 19 ) 

k u(-,t = 0) = po x £ R 2 , 

and the associated free energy takes the form 

•^pksM := / ^logndx / uv dx H — / |x| 2 u dx . (20) 

Jm 2 2 Jb. 2 2 jr 2 
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If (u, v) is a smooth solution of the rescaled Patlak-Keller-Segel system (|19p which decays 
sufficiently at infinity, then 



If we keep in mind the gradient flow interpretation (which is true also for the rescaled equa- 
tion), we can imagine that the limit when t goes to infinity of a solution to the rescaled 
Patlak-Keller-Segel system ()19|) cancel the free energy dissipation ()20p . So that the limit 
solution satisfies the Gelfand equation (fT8]) . 

The question of the speed of convergence is not fully understood. In the case of small mass, 
a linearisation method reduces the problem to a perturbation of the Fokker-Planck equation 
and gives: 

Theorem 7 (Rate of convergence, |12j). Under the assumptions in Theorem^ there exists 
a positive constant M* such that, for any initial data po 6 L 2 (u^ dx) of mass M < M* , 
the rescaled Patlak-Keller-Segel system (fT9|) has a unique solution u G C°(M + , L 1 (M 2 )) n 
L°°((t, oo) x M 2 ) for any r > 0. Moreover, there are two positive constants, C and 5 , such 
that 



As a function of M , 5 is such that limA./->o+ 5{M) = 1 . 

Open question: The question of the speed of convergence for larger mass remains open, 
even if numerical evidence indicates that it should be the case, see [8] . 

2.4. Critical case. In the case M = 8tt, the free energy J-pks is the same as the functional 
which appears in the logarithmic Hardy-Littlewood-Sobolev inequality (jlOj) . The remainder 
entropy which was controlled in (jlip is thus entirely "eaten" by the logarithmic Hardy- 
Littlewood-Sobolev inequality ([TO]) . In [11], we use a three-step procedure: 

How would it blowup: the space is split into balls and annulus. Using the diffusion it 
is possible to prove that in a ball the mass is less than 8tt and to control the influence 
of the interactions outside the ball. So that only when all the mass is concentrated in 
a point, we cannot extend the solution to a bigger interval. If the solution blows up, 
it blows up as a Dirac mass concentrated in the centre of mass. 

When would it blowup: in this case M = 8tt, by the virial computation ([2j), the 2- 
moment remains constant. A De la Vallee-Poussin's type argument, shows that the 
concentration cannot occur in finite time. If the solution blows up, it blows up as a 
Dirac mass concentrated in the centre of mass at infinite time. 

Does it blowup: Keeping in mind the gradient flow structure described in the previous 
section, we can imagine that the solutions tend to the minimisers of the logarithmic 
Hardy-Littlewood-Sobolev inequality (|10p . But as the second-moment is constant 
thanks to the virial computation ([2]), the solution converges to the only minimiser g\ 
of the logarithmic Hardy-Littlewood-Sobolev inequality (|10p which is of finite moment: 
the Dirac mass. 
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As a consequence, we prove 

Theorem 8 (Infinite Time Aggregation, [H]). If the 2-moment is bounded, there exists a 
global in time non-negative free-energy solution of the Patlak-Keller-Segel system ([1]) with 
initial data pq. Moreover if {i p } p eN — oo as p —> oo, then t p \— > p(t p , x) converges to a Dirac 
of mass 8 ir concentrated at the centre of mass of the initial data weakly-* in the sense of 
measure as p — > oo. 

In the radial and bounded case, the blowup rate and refined asymptotics estimates are 
given in the following 

Theorem 9 (Blowup profile, [37] )• In the radial case, in the ball, consider a solution to the 
Patlak-Keller-Segel system ([T]) with dp/dv = pdc/dv and c = on the boundary. Then when 
t goes to oo, 



In [52], C. Sire and P.-H. Chavanis predicted this result by a formal argument considering 
only the first order correction terms. In [37], N. Kavallaris and P. Souplet study the Patlak- 
Keller-Segel system Jl| on a bounded domain with Dirichlet condition, they make successive 
appropriate change of variable and of function, to reduce to a degenerate parabolic problem, 
so that the precise results are difficult to translate back to the original Patlak-Keller-Segel 
system ([1]). In particular, they prove that the solution is the sum of a quasi-stationary profile 
and of a correction term which is significant only for x bounded away from 0. 

The extension of Theorem [8] to the case when the second moment is not finite allows the 

solution to converge to the other minimisers of the logarithmic Hardy-Littlewood-Sobolev 

inequality ()10p . For this purpose we need to introduce another free energy functional which 

still has to be fully understood. Let us first recall the Fokker-Planck version of the fast 

du 

diffusion equation corresponding to the fast diffusion equation — = Ay/u by a self-similar 

at 

change of variable: 



This equation can also be written in a form analogous to ([9]): following [40] for A > 0, define 
the relative entropy of the fast diffusion equation with respect to the stationary solution g\ 






(21) 



by 




2 



Equation (f2Tj) can be rewritten as 
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with 

5H\[u] _ J_ J_ 

But the connection can be seen through the minimisers of T~L\ which are the same as those 
of the logarithmic Hardy-Littlewood-Sobolev inequality (fTU|) . The functional H\ is a weighted 
distance between the solution and q\ and its unique minimizer is q\. It is tempting to compute 
the dissipation of Tl\ along the flow of solutions to the Patlak-Keller-Segel system (fTJ): Let p 
be a sufficiently smooth solution of the Patlak-Keller-Segel system ([1]). Then we compute 

d r , *, 1 /" |Vp| 2 /" o /9 /Mvr f M\ 

In the critical case M = 8-7T the dissipation of the %\ free energy along the flow of the 
Patlak-Keller-Segel system ^ is 

We use the following Gagliardo-Nirenberg-Sobolev inequality due to J. Dolbeault and M. Del 
Pino, see |28j : For all functions / in M 2 with a square integrable distributional gradient V/, 

7T / |/| 6 dx< / (Vjfdx I |/| 4 dx , 

Jr 2 Jr 2 Jr 2 

and there is equality if and only if / is a multiple of a translate of g 1 ^ 4 for some A > 0. 

As a consequence, taking / = p 1 / 4 so that J R2 / 4 (^) dx = 8tt, we obtain T>[p] > 0, and 
moreover, T>[p] =0 if and only p is a translate of q\ for some A > 0. 

Remark 1. While I was finishing this survey, J. A. Carrillo and I realised that this new free 
energy gives another proof of non existence of global-in-time solutions in the super-critical 
case. Indeed, by V[p\ > and (g2J, < H\[p(t)] < 4 V / Mvr/A(l - Af/8vr)t. So that in the 
case M > 8ir, if l-L\[po\ is initially bounded, and thus of infinite 2-moment, then there cannot 
be global-in-time solutions. 

Based on this remark, the main results of [9] in the critical case can be summarised in the 
following: 

Theorem 10 (Existence of global solutions, [9]). Let po be any density in M 2 with mass 8tt, 
such that JtksM < oo, and for some A > 0, W\[p] < oo. Then there exists a global free 
energy solution of the Patlak-Keller-Segel equation (|3U|) with initial data pq. Moreover, 

lim Jpks[p(*)] = FpksIqx] and lim \\p(t) - g x \\i = . 

t— >oo t— >oo 

We even prove further regularity using the propagation of the L p -estimates and the hyper- 
contractivity property of the equation. This theorem can be translated in terms of the Wasser- 
stein distance of the initial data to £>a thanks to the Talagrand inequality: 
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If we keep in mind that the 2-moment can be seen as the Wasserstein distance between the 
solution and the Dirac mass, we see that Theorem [10] completes the picture which emerged 
from the first result set out at the beginning of the section. As soon as we start at a finite 
distance from one of the minimisers g\ we can construct a solution which converges towards it. 
Note that this result is true for the solution that we construct as we do not have uniqueness of 
the solution to the Patlak-Keller-Segel system, even if we strongly believe that this is the case. 
Also observe that the equilibrium solutions g\ are infinitely far apart: let tp(x) = y / A//i|x| 2 /2, 
one has V^?#^ = Q\. Thus, 

Qfj,(x) dx = +00 



2 Jr2 y p 



since the equilibrium densities q\ all have infinite second moments. In particular, "H\[Qu\ = 
+00 for fx 7^ A. There may still be initial data out of these basins of attraction. 

Concerning the proof of Theorem 1 101 we expect the propagation of the bounds on .FpksM 
and T>[p] to give compactness. Unfortunately, T>[p] is a difference of two functionals of p that 
can each be arbitrarily large even when T>[p] is very close to zero. Indeed, for M = 8ir and 
each A > 0, T>[q\] = while 

lim IIoaIU/2 = 00 , lim HVoV 4 ^ = 00 , and lim 5\ = 8ir5o . 

Likewise, an upper bound on J-pks [p] provides no upper bound on the entropy L 2 P log p- 
Indeed, J-pksI/ 3 ] takes its minimum value for p = Q\ for each A > 0, while lim^o J Qx log ^A = 
00. Fortunately, an upper bound on both H\[p] and J^pksI/?] does provide an upper bound 
on f p log p: 

Theorem 11 (Concentration control for J^kSj [S]) - Let p be any density with mass M = 8n 
such that H\[p] < 00 for some A > 0. Then there exist 71 > and Cccf > depending only 
on A and T-L\[p] such that 



Ti 



/ p log p dx < J"pks [p] + Cccf 
Jr 2 



Here we also prove that since T-L\ controls concentration, a uniform bound on both T~L\ and 
T> does indeed provide compactness. We shall prove: 

Theorem 12 (Concentration control for D, [9]). Let p be any density in L 3 / 2 (IR 2 ) with mass 
8tt such that J-~pks[p] is finite, and H\[p] is finite for some A > 0. Then there exist constants 
71 > and Cccf > depending only on X, %\[p\ and J-pksI/ 3 ] such that 

72 f |V/ 4 | 2 dx < irV[p] + Cccd • 

The proofs of this two theorems leads on the following lemma: 

/ yj\ + \x\ 2 p(x) dx < 2\f\M + 2M 3 / 4 (A/^) 1/4 yJU x [p] ■ (23) 
Jr 2 

As explained at the beginning of the section, in [11] we managed to find a ball in which 
the mass was smaller than 8tt. Here, (f23|) gives a vertical cut to prove Theorem [TTJ Indeed, 
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we split the function p in two parts: given /3 > 0, define pp(x) 
(5 large enough, p — pp is such that: 



min{p(a;) , 0}. Bv (|25]> . for 



(p-P/3)<y + C 2 VHlR<^ + 



87T — £q 



< 8 7r — e 



We then apply the logarithmic Hardy-Littlewood-Sobolev inequality method as in (|12p to the 
function p — pp whose mass is less than 8tt. 

The same idea works for the Gagliardo-Nirenberg-Sobolev inequality to prove Theorem 1121 
Let / := /9 1//4 , we split / in two parts by defining fp := min{/, /3 1 / 4 } and hp := / — fp. We 
use (|23p and apply the Gagliardo-Nirenberg-Sobolev inequality to control fyg. 

Idea of the proof of Theorem HOt It follows the line of the convergence of the JKO 
minimising scheme (|14p exposed in the previous section to obtain the Euler-Lagrange equa- 



(24) 



tion ([To]) . Dividing the Euler-Lagrange equation (fT5"j) by y Pt +1 we obtain: 



2V 



,n+l 



Vg 



n+1 



,ra+l 



where V</?"#/9™ +1 = p™. Integrating (|24p we obtain 



n+l Vc n+l 



2Vv> 



2 . 

dx = 




Jm 2 


r 



„rc+l 



dx = W 2 (^,p! 



which is bounded thanks to (|17|) , But the left hand side is a sum of two terms so that 
we cannot conclude any compactness on each of them. That is the reason why we have to 
introduce a regularisation of the Green kernel: let 7 be the standard Gaussian probability 
density in M 2 . Then, for all e > define %(x) = e _2 7 (x/e), and define the monotonically 
regularised Green function G e = 7 e * G * j e . And the corresponding regularised free energy 
functional defined, for all < e < 1, by 



PKS 



[p] '■= [ P(x) \ogp(x) dx- - [[ 
Jr 2 1 J Jm 



p(x) G e (x - y) p(y) dx dy . 



The regularised chemical attractant density c e 
the triangle inequality, from (|24p we obtain 



M 2 xl 2 

G t * p{x) is uniformly bounded, so that by 



zp?srph 



< 



< 



1/2 l 

\Vc e (x)\ p(x) dx ) -| — 
1/2 1 

|Vc £ (x)| 2 p(^) dx J +-W 2 (p,Po)- 



\x — V<p(x)\ 2 p(x) dx 



1/2 



Which gives a bound, which depends on e, on HVy^lb and so on || p ||p for all p £ (1, 00) by 
the Gagliardo-Nirenberg-Sobolev inequality ([5]). The uniform estimates are of course more 
difficult to get. We could probably prove the existence by using the regularisation procedure 
but it is more elegant to prove it following the JKO minimisation scheme (|14p . This optimal 
transport framework is better adapted because in this case we can use displacement convexity. 
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The displacement convexity of H\ is formally obvious from the fact that 

nx[u]= L (~ 2a ^ + v^^ (x) ) dx+c - 



where — yju(x) and |x| 2 u(x) are displacement convex. Actually, we can not expand the terms 
as none of them are finite. Thus, we are forced to introduce a regularisation of T~L X . While 
there are many tools available to regularise functions that are convex in the usual sense, there 
is no general approach to regularising functionals while preserving, or at least not severely 
damaging, their formal displacement convexity properties. The following regularisation is one 
of the cornerstones of this proof: 

T-ixA U \ = / /- , r dx 

JR2 y/Qx + 

which is 7,5-displacement convex with 75 < 0. 

At each step, the main a priori estimate comes from the convexity estimate of the type 

%a[poJ - K\[Pi\ > hmsup , 

which give in this optimal transport framework the above the tangent formulation: 



um +1 ]-uM]>\ [ 



1 V P n T 

X ■ T 



(Vtp(x) - x) p™ dx . 



2A {p nf/2_ 

where Vip is such that VV>#p™ +1 = p™- Inferring the Euler-Lagrange equation (fT5"j): — Vp™ +1 + 
/)" +1 Vc™ +1 = (id — Vtp) p™ +1 /t, we obtain a discrete version of the entropy /entropy dissipa- 
tion inequality 

n x [p n T +1 ]-n x [p n T ]<-rV[p n T ]. (25) 

This inequality is a skeleton version of the crucial estimate which allow to apply the standard 
entropy/entropy dissipation method to study the asymptotics. 

We have to pass to the limit in three parameters: 5, e and r. The <5-regularisation was only 
introduced to justify the integration by parts in (|25p . While passing to the limit in 5 we lose 
slightly in the convexity estimate (|25|) . We obtain that there exist A such that 

H x {p} <n x [po]-rV[p] + TA\h\U /3 and U x [p] < H x [po} ~ r V[p] + t^IA ||p|| 4 / 3 . (26) 

However, a tricky interplay between e and r allows to conserve the estimate at the limit. Let 
us give an idea on how to do it at the first step: define Qo > 0, Tq > by 

Qo := Cp -%a[po] and Tq := mini ®° , lj , (27) 



MMU/3 
and then en which depends on r by 

rV3^[ 87r V3^ 7 -2/3^^ o+To ^ CCD) 2/3] = ^ T 2 2 -. (28) 

By ([26]) , our choice of r and Qo in (f27|) implies that 

n x [p] < n x [p ] - rV[p) + ^ = C po -Q - rV[p] + ^< Cpo - rV[p] . (29) 
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On one hand, the Gagliardo-Nirenberg-Sobolev inequality (|8|) induces V[p] > so that ([29]) 
implies that the energy estimates on p needed for the two concentration controlled inequalities, 
Theorems [TT] and [T2] propagate in the sense that 

^pks[/o] < +oo , H\[p] < C po . 
On the other hand, since 7-L\[p] cannot be negative it implies 

V[p] < 9p± . 
r 

The concentration controlled inequality, Theorem 112] thus gives a bound on ((Vp 1 / 4 ^: 

[ Vp 1 / 4 2 dx < - [nV[p] + Cccd] < -— WC po + t*C CC d] ■ 
Jm 2 72 r 72 

By the Gagliardo-Nirenberg-Sobolev inequality we obtain a bound of the form 



p i/2 dx<8 
And, by Holder's inequality, 



Vp 



1/4 



1 



c 3 



dx < [ttC P0 + t^Cccd] ■= — 

r 72 r 



4 / 3 dx 



vl/3 



p 1 / 3 ^ dx < (8vr 
Now using this bound in (|26|) . we obtain 

Hx\p] ~ n x [ P0 } < -rV[p] + r 1 / 3 ^ [^(Stt) 1 / 3 ^ 273 



2/3 



7T 

MA 



W 2 2 (po,p). 



Our choices of e in (|28|) ensures that we do not lose too much in the convexity estimate: 

n x [p]-n x [p ]<-rV[p] + ^T 2 2- e . 

We see here that the extra term Qqt 2 2~ £ / '4 in the convexity estimate remains bounded in 
the sum. It actually converges to when r goes to 0. These proofs are technical and, even 
if some are interesting from the general theoretical point of view and could be used for other 
problems. Moreover, we can propagate the L p -estimates and prove hyper-contractivity and 
to obtain the convergence of the JKO minimising scheme (j!4[) toward a free energy solution 
to the Patlak-Keller-Segel system ([1]). The proof of the asymptotics result is then a standard 
entropy method as set out at the end of Section 12.31 

Very recently, in [20] . E. Carlen and A. Figalli use a argument of Bianchi-Egnell to obtain a 
quantitative stability for the logarithmic Hardy-Littlewood-Sobolev inequality (]10p and prove: 



\\p(t) - £>a||li(i: 2 ) < 



C 



V lo §(e + 1) 
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3. The Non-linear parabolic-elliptic Patlak-Keller-Segel system 
3.1. The model. In higher dimensions the critical quantity is no longer the mass but the 

L d/2 

-norm, see [27]. We can however replace the linear diffusion with a homogeneous non- 
linear diffusion: 

^(t,x) = dw[Vp m (t,x) - p(t,x)V<p(t,x)} t>0,x£R d , 

- A0(i, x) = p(t, x) , t > , x G R d , ( 30 ) 

p(0,x) =po(x) xeR d , 

where m G (0, 1) and d > 3. In astro-physics, this system models the motion of the mean 
field of many self- gravitating Brownian particles, including a version of the Chandrasekhar 
equation modelling gravitational equilibrium of poly-tropic stars [23]. This system is then 
known as the generalised Smulochowski-Poisson system, in dimension 3 see [H] and [25] in 
dimension 1 and 2. 
Define 

1 

K(x) = c d and c d 



\x\ d ~ 2 (d-2)a d 

where a d := 2 ir d ^ 2 /T(d/2) is the surface area of the sphere S^ 1 in R d . Up to a harmonic 
function <f> = IC*p, so that the system (|30p can be rewritten as a non-local parabolic equation: 

^(t,x) = div[Vp m (t,x) - p{t,x)V(IC* p)(t,x)] t>0,xeR d . (31) 

In this case too the mass is preserved and will be denoted M. 

Let p\(x) := X d p(Xx) with A > 0, the diffusion term scales like \ dm+2 A(p™)(\ x) whereas 
the interaction term scales like A 2d div (p\V(K, * p\)) (Ax). Hence the mass-invariant scaling 
of the diffusion term balances the potential drift in (|3ip if 

m = m d =:2\l - M G (1,2). (32) 

This difference of balance was studied to obtain 

Theorem 13 (First criticality, [M1E5]). Let m d be as defined in (f32|) . 

• if m > m d then the solutions to (|30p exist globally in time, 

• if m < m d then solutions to (|3U|) with sufficiently large initial data blowup in finite 
time, 

• if m = m d there exist two constants M\ > and Mi > M\ such that 

— if M < Mi then the solutions to (|3U|) exist globally in time, 

— if M > M2 there exist initial conditions such that the corresponding solution 
blows up in finite time. 

The proof of this theorem relies on the Gagliardo-Nirenberg-Sobolev inequality and is not 
sharp as we will see in Theorem [TJ] below. In dimension 2, corresponding results were obtained 
by [39} I17j. When no confusion is possible, the index d in m d will be omitted and the critical 
exponent will be denoted m in the sequel of this article. 
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/iRd m — 1 2 

which is related to its time derivative along the flow of f)31 1) by 

\ m — 1 



The analogous of the free energy used in the previous section is: 

t^g[p(t)]:= f P^Mi_ l (( JC(x-y)p(t,x)p(t,y)dxdy 



dx . 

3.2. The sub-critical case. In |10j . the functional inequality used is a variant to the Hardy- 
Littlewood-Sobolev (VHLS) inequality: for all h G L 1 (M d ) nL m (R c! ), there exists an optimal 
constant C* such that 

h.^o { J jR d xR d f - y\ 

We define the critical mass by 



(m - l)C*c d _ 

Theorem 14 (Global-in-time existence, |1CH I58j). If uq is of mass M < M c then there exists 
a global weak solution with initial condition uq. Moreover, this solution satisfies the free 
energy/free energy dissipation inequality. 

The proof of existence follows the lines of the 2d Patlak-Keller-Segel system. Indeed, as a 
direct consequence of the VHLS inequality, for any solution p to the nonlinear Patlak-Keller- 
Segel system ([3T]) 

^ (M^ d - M^ d ) \\p(t)\\™ < G[p(t)} < g[p ] < oo . 

In the case M < M c , it gives the concentration controlled analogous to the entropy a priori 
estimate (|12p of the previous section. It should not be difficult to prove the existence of 
global-in-time solutions using the JKO minimising scheme. 

Open question: The convergence to the self-similar solution is still open. By doing the 
porous medium scaling, we can prove, see [TOj, Theorem 5.2], that for any given mass M < M c 
there exists a unique minimiser Wm to the rescaled free energy. Moreover this minimiser 
is non-negative, radially symmetric and compactly supported. We expect this minimiser to 
attract all the solutions but we have not been able to prove it. 

3.3. The critical case. The balance in the mass-invariant scaling of diffusion and potential 
drift can also be seen in the free energy: If h\(x) := X d h(Xx) then 

g[h x ] = X^ d [ ^1 - X d ^ [ [ | 3_ h{x) %) dx dy 

j R d m i z j j K d xR d \x — y\ 

The diffusion and interaction term balance if m = m^. And in this case 

g[h x ] = x d ~ 2 g[h} . (33) 

This scaling of g has to be compared to the case d = 2, see (|13p . This major difference 
will be a key in the application of the concentration-compactness argument. In the case 
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of the 2d Patlak-Keller-Segel model ([I]), the authors were not able to apply a concentration 
compactness argument due to the rigidity in the scaling of the free energy. Indeed, the scaling 
which preserves the mass also preserves the free energy. Here, as will be shown below, it is 
possible to follow the line of Pierre-Louis Lions's original article. 

The minimisers of Q of mass M c are such that there are R > and z € M. d with 

d/(d-2) 



V(x) 



1 

Rd 



c 



R 



if x e B(z,R), 



(34) 



if x e R d \B(z,R) 

where £ is the unique positive radial classical solution to 

in B(0, 1) with ( = on dB(0, 1) 



AC i m ~ 1 ( 1/(m ~ 1] 



Open question: On the contrary to the 2d Patlak-Keller-Segel system, this proves that 
there exist compactly supported stationary solutions not blowing-up at infinite time. We were 
however not able to prove that they attract some solutions. 

Here, thanks to (|33p we can adapt the concentration-compactness method to prove 

Proposition 1 (How would it blowup, |10j). LetT € (0,oo] and a sequence (t k ) k converging 
to T. If 



lim \\p{t k )\\, 

k— >oo 



OO 



then there are a sub-sequence (tk)j and a sequence (xj)j in M. d such that 



lim 



p{t kp x + Xj ) - — V ( — 



where \ k '■= ||p(ifc)||m m ^ rf ^ and V is the minimiser of Q of the form (|34p with \\V\ 



1. 



The main ingredient of the proof is the following: we set v k {x) := Xf p(tk, Xk x ) so that 
ll^fcllm = 1- By the concentration compactness principle there exists a sub-sequence satisfying 
compactness, vanishing or dichotomy. As already discussed, contrary to the 2d Keller-Segel 
system, here 



lim Q[v k ]= lim \\p(tk)\\- m G[p(t k )}=0. 

k— >oo k — ^oo 



As a consequence 



lim 

k— >oo 



v k (t,x)v k (t,y) 



1 



\x - y 



lim — 

fc-s>oo Cd \ m — 1 



\v k \\Z 



G[vk] > 



Whereas this quantity goes to zero if vanishing or dichotomy should occur. 
Open question: Except in the radially symmetric case, we are not able to say if the blowup 
occurs at the centre of mass or if the blowup escapes at infinity. We could not even rule out 
the possibility that a sub-sequence diverges whereas the other does not. 
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3.4. The super-critical case. The answer was clear for the super-critical case for the 2d 
Patlak-Keller-Segel thanks to the constant sign of the derivative of the 2-moment. Here for 
any p solution to (|30p the virial identity is 

4/ \x\ 2 p(t,x)dx = 2(d-2)g[p(t)] . (35) 
dr jRd 

Theorem 15 (Blowup, [531 E51 EDI [58]). If M > M c , there exist initial data of mass M such 
that the L m -norm of the corresponding solution blows up in finite time. 

In |10j . the proof relies on a procedure directly adapted from [63j: let p be a minimiser of 
the form (j3l"j) and consider 

M 



P0 = M P 



Then, 



AI 



in 



1- 



\p\\ m 



m-l \M Ci 

is negative if M > M c . This result combined with (|35p gives the expected result. 

Open question: We cannot exclude the possibility that solutions with positive free energy 

exist globally in time. 



In |14j a more precise partial answer is proven: the blowup time T being given, we can 
look for solution to (1301) of the form 



1 / x \ 1 
p(t,x) = — — 7* — - and c(t, x) = , $ 



s{t) d \s{t) J v ' ' s{t) d - 2 \s{t) 

where s(t) := [d{T - t)] l / d . 

Theorem 16 (Self-similar blowing-up solutions, [II]). There exists M c S (M c ,oo) such that 
for any M e (M C ,M C ], there exists a self-similar blowing-up solution with a radially symmet- 
ric, compactly supported and non-increasing profile satisfying \\p(t)\\i = ||^||i = M for 
t G [0, T) and ||p(i)||oo — >• 00 as t — > T . 

The method relies on the study of a boundary value problem for the following non-linear 
ordinary differential equation: 

u"(r,a) + - — - u'(r,a) + \u(r, a)^" 1 u(r,a) -1 = 0, r G [0, r max (a)) , 
r 

u(0, a) = a , u'(0, a) = , 

with r max (a) G (0,oo] and p = d/{d — 2). We prove that there are global solutions to this 
problem and that the solutions oscillate around the stationary solution 1, see Figure [TJ The 
solution of Theorem [16] corresponds to that which vanishes and the support corresponds to 
its first zero. The proofs rely on ordinary differential equation tools. For recent results in this 
direction, where asymptotic expansions are performed for a — > 00, see [56] . 
Open question: If a > a c is large enough, u(a, ■) may have several zeros, see Figures 
and each hump corresponds to a solution. It is possible to construct self-similar blowing-up 
solutions of any mass? 
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Figure 1. Behaviour of u(.,a) for a > a c , a = a c and a < a c . 

Open question: The stability of blowing-up solutions is also of interest but seems yet unclear 
according to numerical simulations performed in [5_3j. 




Figure 2. Positivity set of u(.,a) with two (a = 50, left) and three (a = 90, 
right) connected components (d = 3). 

4. Concluding remarks 

The Patlak-Keller-Segel models have attracted much attention these last years. The lit- 
erature is vast and drastically increasing. Much work is also currently carried out on the 
original Patlak-Keller-Segel system with a parabolic equation on the chemo-attractant where 
the authors prove global-in-time existence for a mass less that 8tt |18| but also that there ex- 
ist global-in-time solutions for larger masses [5]. There exist many variants of the presented 
models with prevention of overcrowding [62\ [To] or with non- linear chemo-sensitivity [2J 2JJ , 
etc. The Patlak-Keller-Segel model has recently been used as a basis for more complete mod- 
els |19| . This review is dedicated to the parabolic-elliptic Patlak-Keller-Segel model and tries 
to describe the progress made through energy and functional inequalities methods in the idea 
of [62l E EI] • For more complete reviews see [351 ES EQl [33] • 
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From the author's point of view, the most challenging question is the understanding of the 
blowup. And in this direction progress is still to be made. We are now at a point where we need 
to develop new methods to address those questions. The answer could come from interaction 
with the non-linear Schrodinger equation (NLS) and the unstable thin-film equation (UTF). 
Indeed, the Patlak-Keller-Segel, the NLS and the UTF equation have two levels of criticality. 
The first level is given by the homogeneity of the "attractive" and "repulsive" terms in each 
problem. In our particular case, this refers to the aggregation versus diffusion mechanisms. In 
NLS it is the balance between dispersion and nonlinear attraction. As seen above, the balance 
happens precisely for our chosen exponent m = m^. In the NLS equation this happens for 
the so-called pseudo-conformal non-linearity, see [57] or [221 Chapter 6]. In the UTF equation 
this happens in the so-called marginal case, see [31 0]. In theses three equations, a second level 
of critically occurs when the attractive and repulsive are balanced. In that particular case, 
and for the three models, there exists a critical value M c of the mass which is the maximum 
value of the mass below which the solutions exist globally in time, see \22\ [63l I44j for the 
pseudo-conformal NLS equation and [H [51] for the marginal UTF equation. Note that mass 
refers to the total number of particles for the NLS equation and the L 1 -norm for (|31[) and the 
UTF equation. 

Let us also point out that in all these three problems, the virial method is an elegant way 
to prove that there are solutions which blowup above the critical mass, but it does not give 
any hint on the mechanism of the blowup. However, in jH], for the NLS equation, the result 
goes further and clarifies the blow-up for super-critical masses close to critical. 

The collapsing solutions to the NLS equation have the form of a rescaled ground state 
soliton and the blowup profile Q has the form of a rescaled steady-state solution of the 
Keller-Segel system. The scaling of the leading term in the collapsing solutions has the same 
order \JT — t in both NLS and KS. 
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